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Non-monophyly of Rhacophorus rhodopus , Theloderma 
and Philautus albopunctatus Inferred from 
Mitochondrial 16S rRNA Gene Sequences 


YU Guo-hua'?, RAO Ding-qit'* , YANG Jun-xing’* , ZHANG Ming-wang!? 
(1. Kunming Institute of Zoology, the Chinese Academy of Sciences, Kunming 650223, China; 
2. Graduate School, the Chinese Academy of Sciences, Beijing 100049, China ) 


Abstract; Mitochondrial gene fragments of 16S rRNA gene of four species ( Rhacophorus rhodopus, R . reinwardtii , 
Philautus albopunctatus and P . rhododiscus ) from 11 populations were sequenced in this study. Homologous sequences of 
R . bipunctatus, Theloderma asperum, T .corticale and Buergeria japonica were obtained by screening the GenBank 
database. After excluding all gaps and ambiguous positions, aligned sequences were 500 bp in length with 115 variable 
sites and 92 parsimony-informative sites. Using B . japonica as an outgroup, phylogenetic relationships were analyzed using 
Bayesian inference, maximum parsimony and maximum likelihood methods. Our results indicated that neither of 
R.rhodopus and P . albopunctatus were monophyletic at the species level. The population of R . rhodopus from Hainan Is- 
land was more close to R . bipunctatus than to populations of R . rhodopus from Yunnan Province. Furthermore, the popu- 
lations of R.rhodopus from Yunnan Province can be divided into two main lineages. Theloderma corticale and 
P . rhododiscus were clustered together and J’. asperum was nested in P.albopunctatus. We considered that 
P . albopunctatus Liu and Hu, 1962, was the synonymy of 7. asperum Boulenger, 1886, and suggested removing 
P . rhododiscus from Philautus into the genus Theloderma . 
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Rhacophorus reinwardtii , R.rhodopus , and R . bip- anus (Fei, 1999; Fei et al, 2005). Most previous studies 
unctatus belong to the R .reinwardtii species group be- about these species were mainly focused on karyotype or 
cause of the skin fold on the upside or down side of the Ag-NOR Ce.g., Tan et al, 1989; Rao & Yang, 1996; Li 
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& Hu, 1996), and the information about the taxonomy 
and phylogeny of these species is very limited. Wilkin- 
son et al (2002) found that R . bipunctatus was more 
closely related to R . reinwardtii than to R.rhodopus . 
However, Frost (2006) placed R.rhodopus Liu and 
Hu, 1959, into the synonymy of R . bipunctatus Ahl, 
1927, without any discussions . 

The genus Philautus, as currently understood, is 
clearly heterogeneous in terms of morphology, and possi- 
bly also of life-history (Bossuyt & Dubois, 2001). Phi- 
lautus albopunctatus was placed into the P . albopunctatus 
group along with P.asper and P.tubercudatus by Fei 
(1999). However, Philautus asper has been placed into 
the genus Theloderma as T.asperum ( e.g., Taylor, 
1962; Liem, 1970; Bossuyt & Dubois, 2001; Frost, 
2006 ). According to Channing (1989), Wilkinson et al 
(2002) and Delorme et al (2005), the genus Thelo- 
derma was the sister taxa of Nyctixalus. However, 
Frost et al ( 2006 ) found that T. corticale was more 
close to P . rhododiscus first, and then to Nyctixalus . 
This incongruence is probably due to sampling differ- 
ence as Wilkinson et al ( 2002) did not include 
P . rhododiscus , while Frost et al (2006 ) did not in- 
clude T. asperum . 

Thus, the starting points of our work are to re-exam- 
ine the phylogenetic relationships among R . reinwardtii , 
R.rhodopus and R. bipunctatus, and to obtain a better 
understanding of the placement of 
P . albopunctatus and P . rhododiscus , although both of 
the species were still placed in Philautus by most her- 
petologists (Fei, 1999; Bossuyt & Dubois, 2001; 
Frost, 2006 ). 


phylogenetic 


1 Materials and Methods 


1.1 Samples 

Mitochondrial gene fragments of 16S rRNA se- 
quences from 14 specimens of R . rhodopus , R . reinw- 
ardtii, P .rhododiscus and P . albopunctatus were am- 
plified and 
R bipunctatus, T asperum and T .corticale were ob- 
tained by screening the GenBank database. Although the 
deep relationships among the subfamily Rhacophorinae 


sequenced. Homologous sequences of 


are not clear, the genus Buergeria, distinguished from 
other Asian/African rhacophorids by the absence of a 
bony knob on the third metacarpal ( Wilkinson et al, 
2002 ), definitely is the sister taxon of other Asian/ 
African rhacophorids based on recent molecular studies 
(Richards & Moore, 1998; Wilkinson et al, 2002; Frost 


et al, 2006). Therefore we prudently selected one species 
of Buergeria as an outgroup in the present study. All the 
retrieved sequences from GenBank were: 
japonica, AF458123; R . bipunctatus, DQ283050; 
T corticale, AF458148; T.asperum, AF458144. All 
species used in this study are listed in Tab. 1. 


Buergeria 


Tab. 1 Species used in this work 








GenBank 

Species Sample locality Accession 

Number 

Rhacophorus Simao, Yunnan (SM ) EF646362 

rhodopus Yongde, Yunnan ( YD) EF 646363 

Lvchun, Yunnan ( LC ) EF 646364 

Jingdong, Yunnan (JD ) EF646365 

Xishuangbanna, Yunnan (XSBN) EF 646366 

Longlin, Yunnan ( LL) EF 646367 

Longlin, Yunnan ( LL) EF646368 

Limu Mountains, Hainan (IM) — EF 646369 

Limu Mountains, Hainan (IM) — EF646370 

R. reinwardtii Lvchun, Yunnan (LC ) EF646371 

Lvchun, Yunnan ( LC ) EF646372 

R. bipunctatus = AF458144 

Philautus Jinxiu, Guangxi (JX) EF646373 
rhododiscus 

P. albopunctatus Jinping, Yunnan ( JP ) EF646374 

Jinxiu, Guangxi (JX) EF646375 

Theloderma = AF458148 

asperum 
T. corticale = DQ283050 
Buergeria = AF458123 


japonica 


1.2 DNA extraction, PCR amplification and se- 
quencing 

All voucher specimens used in this study are de- 
posited at the Kunming Institute of Zoology, the Chi- 
nese Academy of Sciences. Total DNA was purified 
from alcohol-preserved muscle or liver using a standard 
proteinase K/SDS digest extraction method followed by 
phenol/chloroform isolation and ethanol precipitation. 
A total of about 550 base pairs corresponding to part of 
the 16S ribosomal RNA gene were amplified via poly- 
merase chain reaction ( PCR ) using the primers 16Sar 
(5'-CGCCTGTITATCAAAAACAT-3' ) and 16Sbr ( 5'- 
CCGGTCTGAACTCAGATCACGT-3' ) ( Palumbi et al, 
1991 ). The lowercase ar and br indicate the amplified 
directions of light and heavy strand, respectively. PCR 
amplifications were conducted in a total volume of 50 
uL using the following cycling conditions: an initial de- 
naturing step at 95°C for 3 min; 40 cycles of denatur- 
ing at 94°C for 1 min, annealing at 51°C for 1 min and 


No. 4 YU Guo-hua et al: Non-monophyly of Three Rhacophoridae Species Inferred from 16$ rRNA Sequences 439 





extending at 72°C for 1 min; and a final extending step 
at 72°C for 10 min. 

PCR products were purified using the “High Pure 
PCR Product Purification Kit” ( Bioteke ) and se- 
quenced in an automated DNA sequencer ( ABI PRISM 
3730 ) using the BigDye terminator v3.1 by Shanghai 
Sangon Biological Engineering Technology & Services 
CO., Ltd. Both sense and antisense strands were se- 
quenced using the corresponding PCR primers. All se- 
quences were deposited in GenBank ( Accession Nos. 
EF646362 - EF646375 ). 

1.3 Molecular and phylogenetic analyses 

Sequences were aligned using CLUSTAL X ver- 
sion 1.83 (Thompson et al, 1997) with default gap 
costs and revised by eye in order to maximize homology 
of position. After cutting off both ragged sides, we ob- 
tained sequences of 540 bp. Considering the uncertain- 
ties of alignment of hypervariabel regions, we preferred 
using a strongly conservative approach, which excluded 
all regions of the gene fragments that could not be 
clearly aligned among all taxa, as recommended by 
Swofford et al (1996 ). Therefore all ambiguous align- 
ments and gaps were removed from the alignment using 
GBLOCKS 0.91b ( Castresana, 2000) with default pa- 
rameters. The revised alignment was 500 bp in length. 

Observed proportional sequence divergence ( p-dis- 
tance ) and the number of transitions and transversions 
in pairwise comparisons were obtained using the com- 
puter program MEGA ( version 3.1; Kumar et al, 2004 ) 
and PAUP* ( version 4.0b10; Swofford, 2002), re- 
spectively. To test for the saturation in base substitu- 
tions, we plotted p-distances against pairwise transition 
and transversion differences. The amount of sequence 
saturation is inferred from the shape of the trend line, 
with a linear relationship indicating that the sequence 
is unsaturated and an asymptotic relationship indicating 
the presence of saturation ( Guo et al, 2005). The sat- 
uration plot was showed in Fig. 1, which confirmed 
that 16S fragments did not reach saturation. 

Prior to phylogenetic analyses, a g1 test was con- 
ducted to assess the amount of phylogenetic signal 
(Hillis and Huelsenbeck, 1992). We generated 10 000 
random trees and calculated the skewness ( g1 ) of the 
resulting tree length distribution with PAUP* 4.0b10. 
The data was subjected to three different phylogenetic 
analyses using: maximum parsimony (MP ), maximum 
likelihood (ML ) and Bayesian inference (BI). MP and 
ML analyses were carried out with PAUP*. The MP 





method was performed using heuristic searches with 
1000 random addition sequence replicates, tree bisec- 
tion reconnection (TBR ), branch swapping and transi- 
tions and transversions given equal weight. Maximum 
likelihood (ML ) analysis was based on the best substi- 
tution model and phylogenetic parameters, which was 
selected using the Akaike Information Criterion ( AIC: 
Akaike, 1973) as implemented in the computer pro- 
gram Modeliest version 3.7 ( Posada & Crandall, 
1998). For the ML analysis, a heuristic search with 
the TBR branch swapping algorithm and 10 random ad- 
dition replicates was employed. The robustness of the 
resulting MP and ML tree topologies was tested using 
bootstrap analyses ( Felsenstein, 1985), with 1 000 
replicates for MP and 100 for ML. Finally, BI was car- 
ried out using MrBayes version 3.1.2 ( Huelsenbeck & 
Ronquist , 2001 ). To ensure that the Bayesian analyses 
were not trapped in local optima ( Huelsenbeck & Boll- 
back, 2001; Leaché & Reeder, 2002), two runs were 
performed simultaneously with four Markov chains 
(three heated and one cold ) starting from random 
trees. The Markov chains were run for 1 000 000 gen- 
erations, sampling every 100 generations thinned the 
data to 10 000 sample points. The program Tracer 1.3 
( Rambaut & Drummond, 2003 ) was used to determine 
when the log likelihood (In L ) of sampled trees reached 
a stationary distribution. Generations sampled before 
the chain reached stationarity were discarded as burn- 
in, and the remaining trees were used to estimate 


Bayesian posterior probabilities ( BPPs ). 


50. 
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2 Results 


2.1 Sequence and tree statistics 

From 10 000 random trees, the g1 statistic signifi- 
cantly supported that our data contained sufficient phy- 
logenetic signals: g1 = -0.96; P <0.01; mean + SD 
tree length = 412.43 + 20.78. Of 500 bp in the align- 
ment, 115 were variable and 92 positions were phylo- 
genetically informative (23.0% and 18.4%, respec- 
tively ). Pairwise uncorrected p distances among taxa 
are shown in Tab. 2. Levels of sequence divergence 
between the outgroup and ingroup lineages ranged from 
9.4% (between B. japonica and R reinwardtii ) to 
12.6% (between B . japonica and P . albopunctatus ). 
Sequence divergence within the ingroup taxa was as high 
as 11.8% ( R reinwardtii compared to P . albopunctatus , 
and Theloderma compared to R.rhodopus ). The 
smallest divergence between two species was 
R . bipunctatus and R .rhodopus from Hainan Island 
(0.4% ). The average transition/transversion ( Ti/Tv ) 
ratio was 3.0 among ingroups. 

The GTR +I + G model was selected as the best-fit 
model of nucleotide substitution for the ML analyses us- 
ing Akaike Information Criterion ( AIC ) as implemented 
in the Modeltest 3.7. Settings for this model were as 
follows: R-matrix = ( 14.5704, 27.8375, 14.2252, 
0.6341, 110.7702 and 1.0000); base frequencies =( A 
= 0.3062, C= 0.2562, G = 0.2196 and T = 0.2180 ); 
proportion of invariant sites = 0.3694; and the shape- 


parameter of the gamma distribution = 0.3840. This 


model was also used to calculate the phylogenetic pa- 
rameters when BI analysis was performed. 

We obtained two most parsimonious trees ( MPT ) 
with 194 evolutionary steps, a consistency index ( CI) 
of 0.768 and a retention index (RI) of 0.855. The 
likelihood values of the ML tree was In L = -1596.63, 
and the likelihood values of the consensus tree in the 
Bayesian approach were In L = — 1609.76 and - 1612.12 
for cold chain of run 1 and run 2, respectively. The 
burn-in in the BI occurred before 250 000 generations, 
so the first 2 500 samples were discarded and the re- 
maining 7 500 sampling trees ( whose log-likelihoods 
converged to stable values ) were used to construct a 


50% consensus tree using the command “sumt” . 
2.2 Phylogenetic relationships 


All four analyses gave similar results, differing on- 
ly where bootstrap support was weak, thus only one phy- 
logeny is presented ( Fig. 2). In all analyses, the ex- 
clusive monophyly of Philautus and Theloderma was not 
supported. Philautus rhododiscus was more closely re- 
lated to T. corticale than to P . albopunctatus with 100, 
78 and 97 for BI, MP and ML support values respective- 
ly, while 7. asperum was nested in P . albopunctatus 
(100, 79 and 91 for BI, MP and ML support values re- 
spectively ). Within R.rhodopus the population from 
Hainan Island was sister to R . bipunctatus with strong 
support (100, 100 and 98 for BI, MP and ML, respec- 
tively ). Populations from Yunnan Province of mainland 
China can be divided into two clades: one was com- 


Tab. 2 Pairwise uncorrected p distances based on partial 16S rRNA gene sequences from eight 


Species used in present study 








Species 1 2 3 4 5 6 7 8 9 10 Wm 12 13 14 «15 «16 © «617 
Rhacophorus reinwardtii 

R. reinwardtii 0.000 

R.rhodopus (LM) 0.060 0.060 

R.rhodopus (LM) 0.060 0.060 0.000 

R. bipunctatus 0.058 0.058 0.004 0.004 

R.rhodopus (YD) 0.066 0.066 0.052 0.052 0.050 

R.rhodopus (JD) 0.066 0.066 0.052 0.052 0.050 0.000 

R.rhodopus (LL) 0.066 0.066 0.052 0.052 0.050 0.000 0.000 

R.rhodopus (LL) 0.066 0.066 0.052 0.052 0.050 0.000 0.000 0.000 

R.rhodopus (LC) 0.076 0.076 0.046 0.046 0.048 0.034 0.034 0.034 0.034 

R.rhodopus (SM) 0.076 0.076 0.046 0.046 0.048 0.034 0.034 0.034 0.034 0.000 

R. rhodopus (XSBN) 0.076 0.076 0.046 0.046 0.048 0.034 0.034 0.034 0.034 0.000 0.000 

Buergeria japonica 0.094 0.094 0.096 0.096 0.098 0.100 0.100 0.100 0.100 0.106 0.106 0.106 

Philautus rhododiscus 0.094 0.094 0.094 0.094 0.096 0.086 0.086 0.086 0.086 0.100 0.100 0.100 0.116 

Theloderma caticale 0.114 0.114 0.108 0.108 0.108 0.110 0.110 0.110 0.110 0.118 0.118 0.118 0.114 0.074 

P. albopunctatus (IX) 0.118 0.118 0.104 0.104 0.106 0.116 0.116 0.116 0.116 0.110 0.110 0.110 0.126 0.096 0.110 

T . asperum 0.116 0.116 0.106 0.106 0.108 0.118 0.118 0.118 0.118 0.112 0.112 0.112 0.124 0.100 0.114 0.006 
P.albopunctatus (P) 0.112 0.112 0.106 0.106 0.108 0.114 0.114 0.114 0.114 0.112 0.112 0.112 0.118 0.100 0.110 0.016 0.014 
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Fig. 2 Bayesian tree based on 500 bp of 16S rRNA with 
B. japonica defined as the outgroup; BI/MP/ML support values 
( % ) are given for nodes that occurred in all applied methods of 
tree reconstruction . 


posed of populations from Southern Yunnan Province in- 
cluding Lvchun, Simao and Xishuangbanna County with 
robust support (100 for BI, MP and ML ), while another 
one consisted of populations from Southwestern Yunnan 
Province including Yongde, Longlin and Jingdong 
County with strong support (100/100/99 BI, MP and 
ML support values respectively ). The average p dis- 
tances within these two clades were 0, and the average 


distance between these two clades was 3.4% . 
3 Discussions 


Wilkinson et al ( 2002 ) found that R. reinwardtii 
was the sister taxon of R. bipunctatus. However, our 
results strongly support that R. reinwardtii is close to 
the clade consisting of R. rhodopus and R. bipunctatus . 
This maybe due to the difference in sequence length. 
Frost ( 2006 ) united R.rhodopus into R. bipunctatus , but 
no discussions were given. According to Fei { 1999 ), 
R. bipunctatus has a green dorsal and two black spots 
on its axilla while the dorsal of R.rhodopus is red- 
brown and has only one black spot on its axilla. Fur- 
thermore, in China, R. bipunctatus is only distributed 
in Xizang while R. rhodopus has a wide distribution in- 
cluding Xizang, Yunnan, Guangxi and Hainan Province . 
Tn our study, the monophyly of R.rhodopus is not sup- 
ported. The population from Hainan Island is the sister- 
group of R. bipunctatus , and the populations from Yunnan 


Province form a clade. Within the Yunnan group, the 
populations from Southem Yunnan with a low elevation 
including counties of Lvchun, Simao, and Xishuangbanna 
form a clade, while the populations from Southwestern 
Yunnan with a high elevation including counties of 
Longlin, Yongde, and Jingdong form another clade. 
Furthermore, the average p distance between these two 
clades is significantly higher than the distances within 
them. These findings indicate that the genetic structure 
of R.rhodopus is complicated. We do not agree with 
Frost ( 2006 ) on placing R.rhodopus Liu and Hu, 
1959, into the synonym of R . bipunctatus Ahl, 1927, 
because the distinctive lineages within R . rhodopus 
may represent some cryptic species and just one of 
these lineages such as the population from Hainan Is- 
land is close to R.bipunctatus. Unfortunately, we 
have no samples from Guangxi Province and adjacent 
countries such as Vietnam and Burma. So more studies 
should be carried to unveil the general phylogenetic 
structure within R . rhodopus . 

The genus Philautus is characterized by the aerial 
direct development of eggs into froglets, without going 
through an aquatic tadpole stage. Although Bossuyt 
and Dubois (2001 ) considered that some of the nomi- 
nal species included in Philautus by Fei (1999 ) do not 
belong to this genus and the whole taxonomy of this 
group proposed by Fei (1999 ) needs re-evaluation, up 
to now there have no demurrers about the taxonomic 
positions of P.albopunctatus and P .rhododiscus . 
Frost et al ( 2006 ) restored P . rhododiscus as the sister 
taxa of T .corticale , but the taxonomy of P . rhododis- 
cus was not changed because, in their study, 
T . corticale and P . rhododiscus were the single repre- 
sentative of Theloderma and Philautus , respectively, 
and the sister relationship between these two species 
only can be treated as an indication of a close relation- 
ship between Philautus and Theloderma . However, in 
this study, the monophyly of Theloderma as restored by 
Wilkinson et al (2002) is not supported. Furthermore, 
P . albopunctatus and P .rhododiscus are not sister taxon 
to each other. As showed in Fig. 2, all analyses supports 
that T .corticale is more closely related to P . rhododiscus 
with high supporting values, and 7’. asperum is nested in 
P . albopunctatus . Moreover, the genetic distance ( uncor- 
rected p distance ) between P . rhododiscus and T. corti- 
cale is as high as 7.4%, and the distance between 
T. asperum and P .albopunctatus is only 0.6%-1.4%. 
Therefore, to correct the paraphyly situation, we suggest 
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that P . albopunctatus Liu and Hu, 1962 should be 
placed into the synonymy species of T. asperum 
Boulenger, 1886, and P . rhododiscus be removed from 
Philautus and put into the genus Theloderma . 


References: 


Akaike H. 1973. Information theory as an extension of the maximum 
likelihood principle [ A ]. In: Petrov BN, Csaki F. Second Interna- 
tional Symposium on Information Theory [ M ]. Akademiai Kiado, 
Budapest. 267-281. 

Bossuyt F, Dubois A. 2001. A review of the frog genus Philautus Gis- 
tel, 1848 ( Amphibia, Anura, Ranidae, Rhacophorinae )[ J]. Zey- 
lanica, 6 (1): 1-112. 

Castresana J. 2000. Selection of conserved blocks from multiple align- 
ments for their use in phylogenetic analysis [ J ]. Molecular Biology 
and Evolution, 12: 814-822. 

Channing A. 1989. A re-evaluation of the phylogeny of old world tree 
frogs [J]. South African Journal of Zoology, 24: 116-131. 

Delorme M, Dubois A, Grosjean S, Ohler A. 2005. Une nouvelle classi- 
fication generique et subgenerique de la tribu des Philautini ( Am- 
phibia, Anura, Ranidae, Rhacophorinae )[ J ]. Bulletin Mensuel de 
la Société Linnéenne de Lyon, 74: 165-171. 

Fei L. 1999. Atlas of Amphibians of China [ M ]. Zhengzhou: Henan 
Publishing House of Sciences and Technology, 1-432. [#% 3%. 
1909. PEI ayy Ae. AGN: WAER AEE. 1- 

432. ] 

Fei L, Ye CY, Jiang JP, Xie F, Huang YZ. 2005. An Illustrated Key to 
Chinese Amphibians [ M ]. Chengdu: Sichuan Publishing House of 
Science and Technology, 1-340. [# 2, MEE, UEF, H 

fE, BOAR. 2005. HEP si RA. WE: WO) BL 
SHAR AEE, 1-340. ] 

Felsenstein JP. 1985. Confidence limits on phylogenies: An approach 
using the bootstrap [J]. Evolution , 39: 783-791. 

Frost DR. 2006. Amphibian Species of the World: An Online Refer- 
ence. Version 4 (17 August 2006 ). The American Museum of Natu- 
ral History, New York, USA. 

Frost DR, Grant T, Faivovich J, Bain RH, Haas A, Haddad CFB, De 
SA RO, Channing A, Wilkinson M, Donnellan SC, Raxworthy CJ, 
Campbell JA, Blotto BL, Moler P, Drewes RC, Nussbaum RA, 
Lynch JD, Green DM, Wheeler WC. 2006. The amphibian tree of 
life[ J]. Bulletion American Museum of Natural History, 297: 1- 
370. 

Guo XG, He SP, Zhang YG. 2005. Phylogeny and biogeography of Chi- 
nese sisorid catfishes re-examined using mitochondrial cytochrome b 
and 16S rRNA gene sequences [J]. Molecular Phylogenetics and 
Evolution , 35: 344-362. 

Hillis DM, Huelsenbeck JP. 1992. Signal, noise, and reliability in 
molecular phylogenetic analyses [J ]. Journal of Heredity, 83: 189- 

195. 

Huelsenbeck JP, Bollback JP. 2001. Empirical and hierarchical Bayesian 
estimation of ancestral states [J ]. Systematic Biology, 50: 351-366. 

Huelsenbeck JP, Ronquist F. 2001. MrBayes: Bayesian inference of 
phylogeny [J]. Bioinformatics , 17: 754-755. 

Kumar $, Tamura K, Nei M. 2004. MEGA3: Integrated software for 


Acknowledgements: We would like to thank LI 
Jia for her help during the laboratorial work. We also 
thank Professor LV Shunging for his tissue samples 
from Hainan Island. 


molecular evolutionary genetics analysis and sequence alignment 
[J]. Briefings in Bioinformatics, 5: 150-163. 

Leache’ AD, Reeder TW. 2002. Molecular systematics of the eastern 
fence lizard ( Sceloporus undulatus ): A comparison of parsimony, 
likelihood, and Bayesian approaches [ J ]. Systematic Biology, 51: 
44-68. 

Li SS, Hu JS. 1996. The karyotypes and Ag-NORS of four batrachian 
species [J ]. Chinese Journal of Zoology, 31(6): 13-16. [ÆR 
YR, HEE. 1996. MEERA ARERR. AE y 
5, 31(6): 13-16. ] 

Liem SS. 1970. The morphology, systematics, and evolution of the old 
world tree frogs ( Rhacophoridae and Hyperoliidae )[ J ]. Fieldiana : 

Zoology , 57: 1-145. 

Palumbi SR, Martin A, Romano $, Owen MacMillan W, Stice L, 
Grabowski G. 1991. The Simple Fool’s Guide to PCR[ M ]. Depart- 
ment of Zoology, University of Hawaii, Honolulu. 

Posada D, Crandall KA. 1998. MODELTEST: Testing the model of 
DNA substitution [J ]. Bioinformatics , 14: 817-818. 

Rao DQ, Yang DT. 1996. Studies on karyotypes of eight species of Rha- 
cophoridae in Yunnan [J]. Sichuan Journal of Zoology, 15 (3): 
105-110. [ EREFF, HAT. 1996. AT 8 PHR ER R 6k 
RCRA. PONSA, 15(3): 105-110. ] 

Rambaut A, Drummond AJ. 2003. Tracer v1.0.1 [CP]. Available from 

http: //evolve.zoo.ox.ac.uk/. 

Richards CM, Moore WS. 1998. A molecular phylogenetic study of the 
old world tree frog family Rhacophoridae [ J ]. Herpetological Jour- 
nal, 8: 41-46. 

Swofford DL. 2002. PAUP* : phylogenetic analysis using parsimony ( * 
and other methods ), version 4.0 [ CP ]. Sinauer Associates. Inc, 
Sunderland, MA, USA. 

Swofford DL, Olsen GJ, Waddell PJ, Hillis DM. 1996. Phylogenetic in- 
ference [ A]. In: Hillis DM, Moritz C, Mable BK. Molecular Sys- 
tematics, second Ed [ M ]. Sinauer Associates. Inc, Sunderland, 
MA, USA. 702-707. 

Tan AM, Wu GF, Zhao EM. 1989. The karyotype of the tree frog Rha- 


cophorus reinwardtii ( Boie ) from Xishuangbanna Daizu Autonomous 








Prefecture, Yunnan Province, China [J]. Chinese Herpetological 
Research, 2(2): 33-36. 

Taylor EH. 1962. The amphibian fauna of Thailand [J]. Univ Kansas 
Sci Bull, 43: 265-599. 

Thompson JD, Gibson TJ, Plewniak F, Jeanmougin J, Higgins DG. 
1997. The CLUSTAL X windows interface: Flexible strategies for 
multiple sequence alignment aided by quality analysis tools [ J ]. 
Nucleic Acids Research , 25: 4876-4882. 

Wilkinson JA, Drewes RC, Tatum OL. 2002. A molecular phylogenetic 
analysis of the family Rhacophoridae with an emphasis on the Asian 
and African genera[ J]. Molecular Phylogenetics and Evolution , 24: 
265-273. 


